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ABSTRACT 


The  USAF  ATHENA  is  a  re-entry  test  vehicle  which  has  four  solid 
propellant  stages.  A  midcourse  correction  type  of  guidance  is  used 
to  compensate  for  boost  dispersions  by  adjusting  the  third  stage 
attitude  and  ignition  time.  The  guidance  computations  are  performed 
in  a  ground  based  digital  computer  using  radar  data  obtained  after 
second  stage  burnout-  Because  this  is  a  re-entry  test  vehicle,  the 
guidance  system  must  cause  the  payload  to  meet  various  re-entry 
constraints  as  well  as  the  usual  impact  point  constraint.  The  com¬ 
plete  derivation  of  the  guidance  equations  and  a  description  of  the 
entire  guidance  loop  are  presented  in  this  report. 
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SECTION  I 

INTRODUCTION 


The  purpose  of  the  USAF  ATHENA  vehicle  is  to  deliver  test  payloads  at 
typical  ballistic  missile  re-entry  conditions  for  the  Advanced  Ballastic  Re- 
Entry  Systems  (ABRES)  program.  The  ATHENA  is  launched  from  Green 
River,  Utah,  and  the  payloads  impact  on  the  White  Sands  Missile  Range 
(WSMR)  in  New  Mexico. 

The  ATHENA  vehicle,  consisting  of  four  solid  propellant  stages  which  are 
burned  to  propellant  depletion,  is  launched  from  a  rail  launcher.  The 
launcher  is  oriented  so  as  to  compensate  for  winds,  based  on  meteorological 
data  obtained  just  prior  to  launch .  The  first  and  second  stages  are  fired  in 
an  upward  direction  and  are  controlled  by  a  combination  of  spin  stabilization 
and  fixed  aerodynamic  fins.  This  type  of  control,  which  causes  the  vehicle 
to  fly  near  zero  angle  of  attack,  has  been  used  on  many  sounding  rockets. 
After  second  stage  burnout  the  vehicle  is  despun  and  a  reaction  jet  attitude 
control  system  orients  the  third  and  fourth  stages  in  a  nominal  attitude 
relative  to  a  prelaunch-erected  gyro  reference  system. 

The  attitude  corrections  and  third  stage  ignition  time  required  to  compensate 
for  dispersions  which  occur  during  first  and  second  stage  operation  are  com¬ 
puted  using  radar  data  acquired  after  second  stage  burnout.  These  computa¬ 
tions  are  performed  in  a  ground  based  digital  computer  at  WSMR  and  the 
necessary  commands  are  then  transmitted  to  the  vehicle.  After  the  vehicle 
has  stabilized  in  the  proper  attitude,  it  is  spun  up  to  a  rate  high  enough  to 
inertially  stabilize  the  third  and  fourth  stages  during  thrusting,  after  which 
the  attitude  control  system  is  separated.  The  third  stage  is  then  ignited  at 
the  computed  time  and  the  fourth  stage  ignites  at  a  predetermined  time  after 
third  stage  burnout. 
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The  purpose  of  this  document  is  to  present  the  derivation  of  the  equations 
used  to  determine  the  required  third  and  fourth  stage  attitude  and  ignition 
time  from  the  radar  data  acquired  after  second  stage  burnout. 
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SECTION  2 


PROBLEM  STATEMENT 


The  ATHENA  guidance  problem  is  similar  to  the  classical  midcourse 
correction  problem;  that  is,  the  position  and  velocity  of  the  vehicle  are 
determined  with  the  use  of  ground  based  radars  and  this  information  is  then 
used  to  determine  the  direction  and  time  of  initiation  of  a  velocity  increment 
that  satisfies  given  constraints.  The  magnitude  of  the  velocity  increment 
cannot  be  varied  in  this  case  since  the  solid  propellant  stages  are  burned  to 
fuel  depletion.  The  direction  of  the  velocity  increment  is  controlled  by 
vehicle  orientation  and  fixed  by  spin-stabilizing  the  third  and  fourth  stages 
after  orientation. 

Position  and  velocity  for  the  guidance  computation  need  be  determined  at 
only  one  instant  of  time,  and  for  simplification  this  is  chosen  as  nominal 
third  stage  ignition  time  for  all  flights.  The  time  of  flight  at  which  actual 
radar  data  becomes  available  varies  from  flight  to  flight,  but  it  is  always 
acquired  after  the  vehicle  has  left  the  sensible  atmosphere.  Hence,  a  form 
of  Kepler's  equations  may  be  used  for  predicting  position  and  velocity  at  the 
predetermined  time  from  position  and  velocity  derived  from  the  radar  data 
at  any  exo-atmospheric  read-out  time  (Reference  1). 

In  addition  to  constraining  the  impact  point,  the  corrections  which  are  com¬ 
puted  by  these  equations  must  constrain  the  vehicle  to  an  angle  of  attack  as 
near  null  as  possible  at  a  predetermined  altitude  in  order  to  satisfy  test 
requirements.  However,  if  the  launch  dispersions  require  corrections  that 
would  reduce  the  re-entry  separation  distance  between  the  payload  and  the 
fourth  stage  below  a  critical  minimum,  it  is  necessary  to  substitute  a 
separation  distance  constraint  for  the  angle  of  attack  constraint.  The  method 
used  to  obtain  the  solution  to  the  guidance  problem  is  discussed  in  the  next 
section. 
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SECTION  3 


METHOD  OF  SOLUTION 


The  analytical  description  of  the  problem  is  derived  from  two  basic  sets  of 
information.  These  are  (1)  the  classical  equations  of  motion  and  (2)  the 
equations  which  analytically  express  the  desired  program  constraints. 

The  equations  of  motion  are  written  in  vector  form  in  an  earth-centered  iner¬ 
tial  coordinate  system  and  are  integrated  to  obtain  the  dynamics  equations. 

The  variables  in  these  equations  are  the  body  attitude  (expressed  in  direction 
cosines),  the  third  stage  ignition  time  (expressed  in  terms  of  the  fourth  stage 
burnout  time),  the  time  and  position  at  the  test  altitude  and  impact,  and  the 
velocity  at  the  test  altitude.  The  dynamics  equations  comprise  a  set  of  nine 
nonlinear  equations  with  fifteen  unknowns. 

Six  constraint  equations  express  the  desired  impact  point  and  either  the  angle 
of  attack  constraint  or  separation  distance  constraint  in  terms  of  the  variables 
in  the  equations  of  motion.  This  results  in  a  total  set  of  fifteen  equations 
with  fifteen  unknowns.  The  simultaneous  solution  of  these  nonlinear  equations 
yields  the  required  attitude  and  ignition  time  of  the  third  stage. 

A  linear  solution  of  these  equations,  made  possible  by  expanding  the  nonlinear 
terms  in  a  Taylor  series  or  equivalent  technique  and  retaining  only  first  order 
terms,  will  not  result  in  a  solution  which  is  sufficiently  accurate  over  the 
required  range  of  inputs.  The  significant  part  of  the  inaccuracy  in  this  solu¬ 
tion  is  a  result  of  the  linearization  of  terms  involving  products  of  two  varia¬ 
bles.  Because  all  of  these  terms  are  functions  of  the  body  attitude,  which  is 
determined  in  the  solution,  it  is  possible  to  recompute  the  "nominal"  values 
in  the  linear  approximations  for  these  terms  after  completing  the  solution 
of  the  linearized  equations.  This  leads  to  an  iteration  procedure  which 
converges  on  a  solution  which  is  sufficiently  close  to  the  solution  obtained 
from  a  complete  trajectory  integration  program  to  meet  all  of  the  ATHENA 
requirements. 
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The  solution  therefore  consists  of  the  following  steps. 


O 


a.  Integrate  the  differential  equations  of  motion  to  obtain 
nonlinear  dynamics  equations. 

b.  Write  algebraic  equations  which  analytically  describe  the 
system  constraints. 

c.  Linearize  these  equations,  where  necessary,  using  a  Taylor 
series  expansion  or  equivalent  technique  and  retain  only 
nominal  and  first  order  terms. 

d.  Solve  the  linearized  equations  for  the  unknown  variables, 
which  include  the  direction  cosines  of  the  body  axis  and  the 
third  stage  ignition  time 

e.  Recompute  the  "nominal"  terms  which  are  functions  of  body 
attitude  in  the  Taylor  series  approximation. 

f.  Repeat  steps  d  and  e  until  the  body  attitude  does  not  change. 

The  derivation  of  these  equations  is  described  in  detail  in  the  remainder  of 
this  report.  The  solution  has  been  programmed  on  a  digital  computer  and  is 
presently  part  of  the  real-time  computer  program  that  is  used  in  the  guidance 
and  command  of  the  ATHENA  vehicle. 
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SECTION  4 


THE  DERIVATION 

4.  1  DYNAMICS  EQUATIONS 

The  dynamics  equations  are  derived  in  vector  form  starting  with  Newton's 
law, 


F  =  mD  (1) 

Figure  1  is  a  geometric  description  of  the  problem. 

The  acceleration  of  the  vehicle  results  from  gravitational  and  thrust  forces, 
therefore 


F  =  m(a  +  £)  (2) 

and 

D  =  a  +  £  (3) 


Integrating  both  sides  of  Equation  (3)  from  t^  to  t  (Figure  2  is  a  description 
of  the  time  notation)  where  t  >  tT  and  rearranging  yields 


(4) 
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Figure  1.  Geometric  Description 


8 


Figure  Z.  Time  Notation  Description 


Integrating  both  sides  of  Equation  (4)  from  tR  to  t  where  t  >  tT  and  re 
arranging  terms  gives 


£  -  °R  +  SrI*  •  *r)  +/'t/  i. dt  dti  +  <*  -  *t>/  T  -  dt 

t /-V  Ia 


''0  0 


‘idtd,i 

K  rR 


(5) 


The  velocity  increment  due  to  the  third  and  fourth  stage  thrusting  is 


The  distance  traveled  during  third  and  fourth  stage  thrusting  due  only  to  the 
vehicle  propulsion  system  is 


(7) 


At  the  time  the  vehicle  reaches  the  test  altitude,  t  =  tT  and  Equation  (5) 
becomes 


£x  =  — R  +  -2r(,t  -  *r>  ♦  £2  +  £S(,t  -  V  *f  T/  *£  d*  dli 

rR  R 


(8) 
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At  the  time  the  vehicle  reaches  the  impact  point,  t  =  tj  and  Equation  (5) 
becomes 

/**I  t  *1 

Dj  =  Dr  +  DR(tj  -  tR)  +  AD  +  AV(tj  -  tr)  + J  1  f  dt  dtj  (9) 

JtR  JtR 

An  equation  which  relates  the  velocity  at  the  test  altitude  to  the  other 
problem  variables  is  required  because  of  the  angle  of  attack  constraint. 

At  the  time  the  vehicle  reaches  the  test  altitude,t  =  t^,  and  Equation  (4) 
becomes 


(10) 


Since  the  air  mass  is  assumed  to  rotate  with  the  earth,  the  required  velocity 
equation  is 

4.2t  =  £t-;:x-2t  <"> 

Substituting  Equation  (10)  into  Equation  (11)  gives 


-gj  Dt  =  Dr  +  AV  +  J  T£  dt  -  w  X  Dr 


(12) 


Equations  (8),  (9)  and  (12)  represent  the  nine  scalar  equations  that  define 
the  dynamics  from  tR  to  tj. 
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4.  2  LINEARIZING  THE  DYNAMICS  EQUATIONS 

In  Equations  (8)  and  (9).  terms  of  the  form  (AV)  (t.)  where  i  =  T,  I,  t  must 
be  linearized  because  they  are  products  of  the  unknown  variables.  The 
linearization  in  this  case  is  performed  in  the  following  manner. 

Let 


AV  =  AVn  +  6AV 


(13) 


and 


*T  = 


tTN  + 


6t„ 


(14) 


therefore 


(AV)(tT)  =  (AVn  +  6AV)(tTN  +  6tT)  (15) 

Expanding  and  neglecting  products  of  perturbations 

(AV)(tT)  «  AV(tTN)  +  AVN(tT)  -  A^N(tXN)  (16) 

Similar  expressions  result  for  tj  and  t  . 

In  Equations  (8),  (9)  and  (12)  the  terms  involving  £  are  linearized  as 
follows. 


f  &  dt */  &  dt  "  .StN^TN  "  *T^ 

*R  *R 


i  t 


m 
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and 


f‘T 

FN  ftl  ,.  ..  ,  f  tTN 

I  &dt  dtj  +  (tT  -  tTN) J  &  dt 

(18) 

lR  ' 

*R  lR 

A  similar  equation  results  for  the  double  integral  of  j  between  t^  and  tj. 

The  linearized  form  of  Equations  (8),  (9).  and  (12)  is  presented  in  the 
linearized  equation  summary. 

Because  the  solid  propellant  engines  are  burned  to  fuel  depletion,  the 
magnitudes  of  AV  and  are  constants  which  can  be  determined  before 
flight.  The  spin  stabilization  of  the  third  and  fourth  stages  orients  the  vehicle 
in  inertial  space  and  as  a  result  AV  and  AD  are  parallel  and  are  completely 
defined  by  their  magnitudes  and  the  direction  cosines  of  the  body  roll  axis.  That 
is,  AD  can  be  expanded  as  a  function  of  the  body  axis  direction  cosines  as 
follows 


AD  =  |AD|(DCx)  i+  |  AD  |(DCy)  j*  +  |_AD|(DCz)k  (19) 

The  vectors  13^  and  can  be  determined  from  the  radar  data  at  the  pre¬ 
selected  time  t^.  The  vector  w  is  the  rotation  rate  of  the  earth  and  the 
terms  involving  £  have  been  linearized  as  shown  in  Equations  17  and  18. 

Only  the  three  body  axis  direction  cosines  are  unknown  in  the  vectors  AV 
and  AD.  The  other  unknowns  in  Equations  (8),  (9),  and  (12)  are:  DT,  three; 
t,p,  one;  t  ,  one;  Dj,  three;  tj,  one;  and-g^D^,,  three.  The  final  count  finds 
nine  equations  and  fifteen  unknowns.  Therefore,  six  constraint  equations 
must  be  written  to  complete  the  analytical  description  of  the  problem. 
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DIRECTION  COSINE  CONSTRAINT 


4.  3 


The  first  of  the  six  constraints  relates  the  three  body  axis 
of  which  only  two  are  independent.  The  equation  is 

direction  cosines 

(DCx)2  +  (DCy)2  +  (DCz)2  =  1 

(20) 

To  linearize  Equation  (20),  let 

DCx  =  DC*  ,  +  6DCx 

(21) 

DCy  =  DCy^  +  &DCy 

(22) 

DCz  =  DCz.,  +  6dCz 

N 

(23) 

Substituting  Equations  (21),  (22),  and  (23)  into  Equation  (20),  expanding,  and 
neglecting  products  of  perturbations,  yields 

(DCx)(DCxN)  +  (DCy)(DCyN)  +  (DCz)(DCzN)  a  1  (24) 


This  is  the  linearized  direction  cosine  constraint. 

4.  4  IMPACT  POINT  CONSTRAINT 

The  impact  point  constraint  can  be  written  from  inspection  of  Figure  1.  It  is 
simply 


°I=£l 


(25) 


Because  the  transformation  between  the  earth  fixed  coordinate  system  and 
the  inertial  coordinate  system  is  a  function  of  the  unknown  time  tJ(  the  exact 
values  of  thecomponents  of  Ej  are  not  known  before  flight  in  the  inertial 
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coordinate  system.  To  overcome  this  problem,  the  right  hand  side  of 
Equation  (25)  is  replaced  by  the  first  two  terms  of  its  Taylor  series  expansion 

—I*5  J-iN  +  — IN(tI  "  W  *26) 


in  addition 


—IN  =  Tt  ^IN  +  ^  X  —IN 


(27) 


Because  the  rate  of  change  of  JEjN  in  the  earth  fixed  system  is  zero,  Equation 
(27)  becomes 


H  *— IN 


(28) 


therefore 


(29) 

Equation  (29)  is  the  linear  form  of  the  impact  point  constraint. 

4.  5  TEST  ALTITUDE  CONSTRAINT 

The  test  altitude  constraint  defines  the  altitude  at  which  the  angle  of  attack  or 
separation  distance  constraint  applies.  By  inspection  this  constraint  can  be 
written  as  follows. 


TN 


(30) 


Equation  (30)  must  be  linearized  to  allow  a  linear  solution  of  the  whole  system 
of  equations.  This  linearization  is  effected  as  follows. 
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First,  rewrite  the  constraint  as 


which  can  be  written  as 


Let 


— T  '  --TN  +  — T 


(31) 


(32) 


(33) 


Using  Equation  (33)  and  neglecting  products  of  perturbations,  the  left  hand 
side  of  Equation  (32)  becomes 


D_ 
— T 


n 

— T 


~  D 
— TN 


— TN 


2  — TN 


60, 


(34) 


However 


— TN 


(35) 


therefore 


— T 


— T  85  2  — TN 


— T 


-  D, 


TN 


— TN 


(36) 


Substituting  Equation  (36)  into  Equation  (32)  and  rearranging  the  terms  yields 
the  required  linear  constraint. 


— TN  — TN  ~  —  TN  ‘  — T 


(37) 
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ANGLE  OF  ATTACK  CONSTRAINT 


4.  6 

From  a  re-entry  test  viewpoint,  it  is  desirable  to  have  the  capability  to 
provide  near  null  angle  of  attack  at  a  specified  altitude.  This  requirement 
can  be  fulfilled  by  the  ATHENA  vehicle  without  an  active  attitude  control 
system  on  the  payload.  Due  to  its  simplicity  the  vehicle  cannot  always 
completely  null  the  angle  of  attack  even  with  perfect  system  operation. 
However,  the  small  variations  about  null  that  do  occur  do  not  compromise 
the  test  objectives. 

Because  the  vehicle  is  spin  stabilized  the  body  axis  and  AV  (the  velocity 
added  by  the  third  and  fourth  stages)  are  parallel.  Therefore  if  AV  and 
^  D,p  are  parallel,the  angle  of  attack  will  be  nulled.  From  Figure  3  it 
can  be  seen  that  the  integral  of  £  adds  a  degree  of  freedom  in  the  vertical 
plane  that  is  not  present  in  the  horizontal  plane.  Because  there  is  insufficient 
freedom  in  the  horizontal  plane  the  impact  point  and  angle  of  attack  constraints 
can  conflict.  This  conflict  is  best  visualized  as  follows.  It  can  be  seen  in 
Figure  3  that  the  velocity  increment,  AV,  will  have  to  be  parallel  to  the 
vector  sum  of  and  uj  X  D^,  to  null  angle  of  attack  in  the  horizontal  plane. 

In  most  cases  this  attitude  would  not  result  in  a  proper  impact  Therefore 
the  angle  of  attack  in  this  plane  cannot  always  be  nulled.  This  problem  is 
not  present  in  the  vertical  plane  because  of  the  extra  degree  of  freedom 
resulting  from  the  integral  of  £. 

This  restriction  does  not  compromise  the  usefulness  of  the  system  because 
the  deviations  from  zero  angle  of  attack  in  the  horizontal  plane  are  within 
acceptable  limits.  However,  the  angle  of  attack  constraint  equation  must 
be  written  in  a  manner  that  will  not  result  in  unrealizable  solutions. 

Therefore,  it  is  necessary  to  write  a  constraint  equation  that  causes  only 
the  angle  of  attack  in  the  vertical  plane  to  be  nulled.  The  derivation  of 
this  constraint  equation  follows. 
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HORIZONTAL 


VERTICAL  PROJECTION 


Figure  3.  Velocity  Vector  Diagrams 


The  angle  of  attack  constraint  equation  is  written  from  inspection  of 
Figure  4.  The  vector  is  normal  to  the  vertical  plane  defined  by 
Dtn  and  ^_VN:  that  is, 


-TN  "  — TN  X  — N 


(38) 


Therefore,  4  D^,  X  f  lies  in  this  vertical  plane  and  is  normal 
,  at  —  I  —  1  in 

to  7!t  — T' 

if 


(4-£txatn)  »  (”> 

then  AV  must  lie  in  the  plane  shown  in  Figure  4.  Therefore,  the  projections 
of  Dt  and  AV  in  the  vertical  plane  must  be  parallel.  This  satisfies  the 
angle  of  attack  constraint. 

This  constraint  is  linearized  as  follows.  Write  Equation  (39)  in  an  alternate 
form. 

'tn-  = 0  <40> 


Let 


AV  =  AVn  +  6AV 


(41) 


and 


D„  =  -£•  D_..  +  6  4-  Dn 


d  n  _  d 
clt  — T  ~  Ht  — TN 


Ht^T 


(42) 
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Substituting  Equations  (41)  and  (42)  into  Equation  (40),  expanding  and 
neglecting  products  of  perturbations  yields 

■3F— T  '  (-TNX—  n'  +—  '(clt—  TN  X  ^Tn)*  ^N  '  ('3t  ~TN  X  -Tn)  (43) 

Equation  (43)  is  the  linear  form  of  the  angle  of  attack  constraint. 

4.7  SEPARATION  DISTANCE  CONSTRAINT 

Some  of  the  experiments  require  a  minimum  separation  distance  between 
the  payload  and  the  expended  fourth  stage  at  a  given  altitude  Since  the 
separation  velocity  increment  due  to  the  fourth  stage  retro  rockets  is 
determined  before  the  flight  and  separation  occurs  at  a  fixed  time  after 
burnout,  the  separation  distance  can  be  maintained  by  controlling  the  time 
between  burnout  and  arrival  at  the  test  altitude. 

That  is, 

Kt  =  tT  -  tr  (44) 

This  is  the  separation  distance  constraint. 
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The  test  altitude  constraint  is 


o 


— TN  '  — TN  ~  — TN  ■  — T 


(50) 


The  angle  of  attack  constraint  is 


Tt-7 


— TN  X  — N*  + 


AV 


(3t  — TN  X  — TN^  =  AV 


•N 


(cft  — TN  X  -TN* 


(51) 


The  separation  distance  constraint  is 


K 


T  " 


-  t 

T 


(52) 


L) 
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SECTION  6 


THE  ITERATION  LOOP,  LINEARIZATIONS,  AND  ACCURACY 


There  are  three  types  of  terms  that  have  been  linearized.  These  are,  (1) 
terms  which  are  functions  of  body  attitude  and  time,  (2)  vectors  from  the 
center  of  the  earth  to  the  vehicle,  and  (3)  terms  which  are  functions  of 
gravity. 

The  simple  solution  of  the  set  of  linearized  equations  will  not  provide  a 
sufficiently  accurate  solution  for  the  body  attitude  and  third  stage  ignition 
time  if  the  position  and  velocity  at  t^  are  dispersed  from  the  nominal  values 
by  reasonable  amounts.  This  is  mainly  due  to  the  linearization  of  the  body 
attitude  dependent  terms  and  is  overcome  with  the  use  of  the  following 
iteration  procedure. 

Because  the  body  attitude  (expressed  as  direction  cosiies)  is  determined  in 
the  solution  of  the  equations,  it  is  possible  to  recompute  all  of  the  so  called 
"nominal"  terms  which  are  functions  of  body  attitude  after  the  solution  is 
complete.  Starting  with  values  obtained  from  the  nominal  trajectory  and 
recomputing  all  of  the  "nominal"  terms  that  are  functions  of  body  attitude 
after  the  solution  of  the  equations  leads  to  an  iterative  procedure.  This 
iteration  procedure  converges  on  solutions  which  do  not  contain  errors 
resulting  from  the  body  attitude  dependent  terms.  Perhaps  this  is  more 
easily  understood  by  examining  a  particular  example.  Consider  the  term 
(AV)  (tj).  The  linearization  used  for  this  term  is 

(AV)(tT)  =  AV(tTN)  +  ^YN(tT)  -  AVN(tTN)  (53) 

Because  the  criterion  for  convergence  is  that  the  direction  cosines  used  to 
compute  the  components  of  AV^j  are  identical  to  those  found  from  the  solution 
of  the  equations,  the  components  of  AV  and  AV^  are  equal  when  the  iteration 
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is  completed.  Therefore,  the  first  and  last  terms  on  the  right  hand  side  of 
Equation  (53)  cancel  and,  since  AV  equals  AV^,  the  iteration  process  has 
nullified  errors  due  to  this  linearization.  Similar  reasoning  applies  to  all  of 
the  terms  which  involve  body  attitude.  A  diagram  of  the  iteration  loop  is 
shown  in  Figure  5. 

Terms  containing  integrals  of  £  have  been  linearized  by  expanding  them 
in  a  Taylor  series  and  retaining  only  first  order  terms,  with  the  integrals 
of  £  considered  to  be  only  functions  of  time.  Consequently,  the  first  order 
position  dispersions  were  disregarded  as  were  the  higher  order  and 
cross  product  terms.  The  adequacy  of  this  assumption  has  been  verified  by 
comparing  solutions  from  the  guidance  equations  with  equivalent  solutions 
obtained  from  a  trajectory  integration  program.  It  was  found  that  the  gravity 
integral  approximations  introduce  the  most  significant  errors  into  the  guidance 
equations,  but  these  errors  are  an  order  of  magnitude  smaller  than  those 
introduced  by  other  links  in  the  guidance  loop  (e.g.,  radar  and  attitude  con¬ 
troller  errors)  and  consequently  they  can  be  tolerated. 

Because  of  the  nature  of  the  vectors  from  the  center  of  the  earth  to  the 
vehicle,  the  linearization  of  these  vectors  does  not  significantly  affect  the 
accuracy  of  the  solution.  The  solution  of  the  present  equations  with  the 
iteration  on  the  direction  cosines  results  in  an  accuracy  which  easily  meets 
all  of  the  program  objectives  for  the  ATHENA  missions. 
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Figure  5.  System  Block  Diagram 
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APPENDIX  1 

SCALAR  EQUATIONS 


The  expansion  of  the  vector  equations  in  the  linearized  equation  summary 
to  obtain  the  required  scalar  equations  is  presented  in  this  appendix. 

The  scalar  components  of  Equation  (45)  are  obtained  using  the  expansion 
of  AD  shown  in  Equation  (19)  and  a  similar  expansion  for  AV  and  AV^. 

The  scalar  components  are  shown  in  Equation  (54)  for  i  =  x,  y,  z. 


°Ti  =  °Ri  +  6Ri(tT  *  lR>  +  l^lDCi  +  |AV|<tTN  •  W001 
+  |AV|(I)CiN)(txN  -  tTN)  +  |AV|(DCiN)(tT  -  tr) 

*  Jf™  J[‘V d*.  *  -  w/'1  V 

lR  lR  R 


(54) 


Equation  (54)  can  be  rewritten  in  the  form 


A.  -  (-1)Dt.  +  (C.)tT  +  (F.)DCi  +  (G.)tx 


(55) 


where 


Ai  =  ‘tn 


jf‘TNM>  -  /™  f'l p  «,  -  DRi  +  i>Ri«R 

lD  An  ■'t- 


R  lR  "R 

-  !^v|(DCiN)(tTN-tTN) 


(56) 


A- 1 


ci  =  !^V  |{DCiN)  +  Dr.  +  /tTNgidt  (57) 

lR 

Fi=  l£D|  +  I^LI^TN’W  (58) 

Gi  =  -|AV|(DCiN)  (59) 

The  scalar  components  of  Equation  (46)  are  obtained  in  a  manner  similar  to 
the  method  used  to  derive  Equation  (54)  and  can  be  written  in  the  following 
form,  for  i  =  x,  y,  z. 

a.  -  (-l)Dn  +  (e.)tj  +  (f.)DCi  +  <h.)tT  (60) 

where 

■  ‘IN jf ' V  •  jfIN  jfV  *,  -  *  V* 

lR  lR  R 

-|AV|(DCiN)(tTN-tIN)  (61) 

ei=  !^H°CiN  +  DRi  +  JINg.dt  (62) 

lR 

f.  =  |AD|  +  !AV|(tIN-tTN)  (63) 

h.  =  -|AV|DCiN  (64) 
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The  scalar  components  of  Equation  (47)  are  shown  in  Equations  (65),  (66), 
and  (67). 


it  DTx  -  6Rx  +  l*HDC*  +  / T  -  STNxI'TN  -  *T> 

»  T-* 


-  w  D™,  +  w  D_ 

y  Tz  z  Ty 


3TD Ty  =  6Ry  *  l*ilDCV  +  /tTNeyd'  -  STNyf'TN  ‘  V 

»  n 


+  W  Drp  -  u>  D  — 
x  Tz  z  Tx 


Tt  DTx  =  °R*  +  l"lDC2  *  jf  ™*,  dt  -  8tN,<* 

^  T"» 


TN  "  tT) 


w  D_  +  w  D„ 
x  Ty  y  Tx 


These  equations  can  be  written  in  the  following  form  for  i  =  x,  y,  z. 


(65) 


(66) 


(67) 


ki  =  (" 1}  37 °Ti  +  <ni)DCi  +  (ri)‘T  +  (wi)DTx  +  (si)DTy  +  (vi)DTz  (68) 

where 


=  "°Ri  ‘  j  gidt  +  gTNitTN 

R 

(69) 

=  \*L\ 

(70) 

=  gTNi 

(71) 

A-  3 


(72) 


w  =0 
x 


W  =  -w 

y  z 


w  =  +oo 
Z  y 


S  =  +00 
X  z 


s  =  0 

y 


s  =  -oo 
Z  X 


(73) 


V  =  -CjJ 

x  y 


V  =  +oo 

y  x 


v  =  0 
z 


(74) 


The  scalar  components  of  Equation  (49)  are  shown  in  Equations  (75),  (76), 
and  (77). 


DIx  =  EINx  +  (“yEINz  ‘  “zEINy)(tI  “  W 

DIy  =  EINy  +  (wzEINx  "  WxEINz)(tI  ’  W 

Dt  =  Etm  +  (oo  Et.t  -  oo  ET  .  )(t  -  tTKT) 

Iz  INz  x  INy  y  INx  I  IN 

These  equations  can  be  written  as  follows  for  i  =  x,  y,  z. 

Qi  -  <-J>DIi  + 


(75) 

(76) 

(77) 


(78) 


where 

a  =  -E--.  +  (oo  Etxt  -  w  Et1.t  )tT»T 

x  INx  y  INz  z  INy  IN 

ay  =  _EINy  +  (WzEINx  "  WxEINz)<:IN 

az  =  'EINz  +  (UxEINy  '  WyEINxHIN 

^x  ”  ^WyEINz  ~  WzEINy^ 


(79) 

(80) 

(81) 

(82) 
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v  =  (w  ETVr  -  w  ETVr  ) 
'y  z  INx  x  INz 


v  =  (w  Etm.  -  «  E  M  ) 
lz  '  x  INy  y  INx 


(83) 

(84) 


The  scalar  expansion  of  Equation  (50)  is  shown  below 


— TN  ’  — TN  "  — TN  ’  — T 


(85) 


therefore 


I— TnI  =  (DTNx)DTx  +  (DTNy)DTy  +  <DTNz'DTz 


(86) 


The  expansion  of  Equation  (51)  has  the  following  form 


»VaTDT*  *  <Nyl^DTy  +  +  <T*>DC* 

+  (Ty)DCy  +  (Tz)DCz  =  W 


(87) 


where 


=  <TNy!— 1 

DCz^ 

"  fTNz 

|  AV| 

lDCyN 

(88) 

=  JTNz^  — 1 

IDCxn 

"  iTNx 

|AV| 

IE>Czn 

(89) 

=  iTNxl— 1 

lDCyN 

"  <TNy 

\2L\ 

DCxn 

(90) 

=  <|av|>(«XNi£d 

TNy  *  1 

TNy3tDTNz) 

(91) 

■  <  I5Y.  l)(«TNx 

TNz  ‘  1 

TNz 

3t°  TNx) 

(92) 
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and 


Tz  ~  (  I—  I^TNyHt  °TNx  "  iTNxHtDTNy) 

W  =  jAV|DCxN(iTNz^DTNy  -  iXNy3t  dTNz) 

+  lAV  iDCyN(iTNxgt  dtNz  _fTNzatDTNx) 
+  i^V  l°CzN(iTNy^ dTNx  -  ^XNx3tDXNy) 


*TNx  =  DTNy!^!DCzN 


*TNy  "  DTNz  I  — IDCxN 


'tNz^TNx^I0^ 


N 


-  D 

-  D 

-  D 


TNz 


TNx 


1  AV  |DCy 
[AY  |DCz 


XNyl^!DCx 


N 

N 

N 


(93) 


(94) 


(95) 

(96) 

(97) 


The  scalar,  linearized  equations  are  presented  in  matrix  form  in  Equation 
(98)  on  the  following  page.  Any  symbols  that  are  used  in  this  set  of  equations 
that  are  not  defined  in  the  Glossary  are  defined  by  equations  in  this  appendix. 


If  the  earth  centered  inertial  coordinate  system  is  oriented  with  the  z  axis 

aliened  with  the  earth's  spin  vector  w  =  s  =  v  =  v  =0. 

z  z  x  y 

If  the  separation  distance  constraint  is  required,  the  14th  row  of  the  matrix 
is  replaced  by  the  equation 

Kt  =  tT  -  tT  (99) 
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APPENDIX  2 


GLOSSARY 

The  symbols  used  in  this  document  are  shown  graphically  in  Figures  1  and  2 
in  addition  to  being  defined  in  this  glossary. 

A  line  under  a  letter  (x)  is  used  to  denote  a  vector  quantity. 

A  dot  above  a  vector  (x)  is  used  to  denote  the  time  derivative  of  the  vector 
relative  to  an  inertial  coordinate  frame. 

Two  dots  above  a  vector  (3 1)  are  used  to  denote  the  second  time  derivative  of 
the  vector  relative  to  an  inertial  coordinate  frame. 

The  symbol  x)  is  used  to  denote  the  time  derivative  of  the  vector  x 
relative  to  a  rotating  coordinate  frame. 

The  first  subscript  on  a  vector  is  used  to  denote  the  time  at  which  the  vector 
is  expressed. 

The  subscript  N  is  used  to  denote  nominal  parameters.  Those  nominal 
parameters  which  are  functions  of  body  attitude  are  recomputed  for  each 
iteration  (see  Page  25). 

a  =  the  thrust  acceleration  of  the  vehicle. 

D  =  the  vector  from  the  center  of  the  earth  to  the  vehicle. 

DCx 

DCy  =  the  three  directions  cosines  of  the  vehicle  body  roll  axis  in 

DCz  the  earth  centered  inertial  coordinate  system. 

E  =  the  vector  from  the  center  of  the  earth  to  the  desired 

~  vacuum  inpact  point. 

F  =  total  applied  force. 

c  =  the  acceleration  due  to  gravity. 
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orthogonal  unit  vectors  defining  the  earth  centered  inertial 
coordinate  system. 

convergence  constant. 

the  time  between  fourth  stage  burnout  and  arrival  at  the  test 
altitude  that  satisfies  the  separation  distance  constraint. 

the  vector  normal  to  the  nominal  vertical  plane. 

mas  s . 

the  vector  from  the  desired  vacuum  impact  point  to  the  center 
of  the  radar  coordinate  system. 

the  vector  from  the  radar  coordinate  system  center  to  the 
vehicle. 

time  (see  Figure  2  for  specific  definitions), 
the  perturbation  from  the  nominal. 

the  distance  travelled  during  the  third  and  fourth  stage  thrust¬ 
ing  due  only  to  propulsion  forces. 

the  velocity  increment  due  to  the  third  and  fourth  stage 
thrusting. 

the  time  between  third  stage  ignition  and  fourth  stage 
burnout. 

the  rotation  rate  of  the  earth, 
time  at  test  altitude. 

end  point  for  Keplerian  transfer  of  radar  data  (chosen  as 
nominal  third  stage  ignition  time). 

third  stage  ignition  time. 

fourth  stage  burnout  time. 

time  of  impact. 

time  when  adequate  radar  measurements  are  available- 
second  stage  burnout  time, 
launch  time. 
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